clc;
clear;
A=sym(magic(4));
B=pinv(A),
A*B,
B*A,
norm(A*B*A-A)